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Abstract 

Decay to asymptotic steady state in one-dimensional logistic-like mappings is characterized by considering a phenomenological 
description supported by numerical simulations and confirmed by a theoretical description. As the control parameter is varied 
bifurcations in the fixed points appear. We verified at the bifurcation point in both; the transcritical, pitchfork and period-doubling 
bifurcations, that the decay for the stationary point is characterized via a homogeneous function with three critical exponents 
depending on the nonlinearity of the mapping. Near the bifurcation the decay to the fixed point is exponential with a relaxation 
time given by a power law whose slope is independent of the nonlinearity. The formalism is general and can be extended to other 
dissipative mappings. 
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1. Introduction 


Mappings are often used to characterize the evolution of dy¬ 
namical systems by using the so called discrete time ill]. The 
interest in the subject was increased after the investigation of 
May 12] with direct application to biology After that a 
large number of applications involving mappings were consid¬ 
ered particularly related to physics iSllll , chemistry, biology, 
eng ineering, mathematics and many others 
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The type of dynamics certainly depends on the control pa¬ 
rameter. As it is varied, bifurcations appear changing the dy¬ 
namics of the steady state lHI] , and for specific rates eventually 
leads to chaos |[l9l |2^ . Collisions of stable and unstable man¬ 
ifolds yield the destruction of chaotic attractors iBIIil]. Many 
of these dynamical properties are already known and are taken 
indeed at the asymptotic state. The way the system goes to equi¬ 
librium is generally disregarded just by considering the evolu¬ 
tion for a large transient. 

Our main goal in this Letter is to apply a scaling formal¬ 
ism to explore the evolution towards the equilibrium near three 
types of bifurcations in a logistic-like mapping: (a) transcrit¬ 
ical; (b) pitchfork and (c) period-doubling. Indeed at the bi¬ 
furcation point the orbit relaxes to equilibrium in a way de¬ 
scribed by a homogeneous function with well defined critical 
exponents [23l [2^. Such exponents are not universal and 
depend mostly on the nonlinearity of the mapping and on the 
type of bifurcation. Near a bifurcation, the relaxation to the 
equilibrium is exponential, with a relaxation time characterized 
by a power law 1 I 22 I 1 . Here, two different procedures are used to 
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obtain the exponents. The first one is mostly phenomenologi¬ 
cal with scaling hypotheses ending up with a scaling law of the 
three critical exponents. The second one considers transform¬ 
ing the difference equation into a differential one and solving 
it with the convenient initial conditions. Our analytical results 
confirm remarkably well the numerical data obtained via com¬ 
puter simulation. 

This piece of work is organized as follows. First the map¬ 
ping, the equilibrium conditions and a phenomenological ap¬ 
proach leading to the scaling law is described. Then we discuss 
the critical exponents by transforming the difference equation 
into a differential equation. Moving on we present discussions 
and extensions to other bifurcation when finally our conclusions 
are drawn. 

2. The mapping and phenomenological properties of the 

steady state 

The mapping we consider is written as 

x„+i = Rx„(l - xl), (1) 

where y > 1 , is a control parameter and v is a dynamical 
variable. A typical orbit diagram is shown in Fig. [T]for: (a) 
7=1 (logistic map) and; (b) 7 = 2 (cubic map). 

The fixed points are obtained by solving x„+i = = x* 

and two cases must be considered: (i) 7 is an even number 
or; (ii) 7 is any other value (odd, irrational etc). For case (i) 
there are three fixed points. One is x* = 0, which is stable 
(asymptotically stable) for R e [0,1) and the two others are 
X* = +[1 - 1/7?]^/^, which are stable for R e (1, (2 y)/y). 

The bifurcation at R = 1 is called pitchfork ||25ll2^. For case 
(ii) there are only two fixed points. One is x* = 0, stable for 
R e [0,1) and the other is x* = [1 - 1//?]^/^, being stable for 
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Figure 1: Orbit diagrams obtained for Eq. m considering: (a) 7 = 1 (traditional 
logistic map) and; (b) 7 = 2 (cubic map). Some bifurcations are indicated in 
the figures. 

R € (1,(2 + y)/y). Transcritical is the bifurcation aXR = 1 for 
this case. Both bifurcations are identified in Fig. [TJa^b). Our 
goal is to consider the convergence to the fixed point v* = 0 
at the bifurcation in Rc = 1 and in its neighbouring such that 
/d = Rc - R = 0, with R < Rc. 

The orbit diagram allows also to extract more properties. Af¬ 
ter a transcritical bifurcation, as seen in Fig. Win), the period-1 
orbit is stable for the range R e (1,3) when a period-doubling 
bifurcation happens. After that a period-doubling sequence is 
observed, obeying a Feigenbaum scaling |tl9l |20|] until reach 
the chaos. Similar dynamics, for a different range of R is also 
observed after a pitchfork bifurcation, as shown in Fig. [Hb). 

The natural variable to describe the convergence to the steady 
state is the distance from the fixed point ll23ll . Indeed for the 
fixed point x* = 0 , the distance to the fixed point is the own 
dynamical variable x. The convergence to the steady state must 
also depends on the number of iterations n, on the initial condi¬ 
tion xo, and of course on the parameter /i = Rc-R. The param¬ 
eter fi = 0 defines the bifurcation point and the convergence to 
the fixed point is shown in Fig. [2] for two different values of y: 
(si) y = 1 and; (b) 7 = 2 and different initial conditions xq, as 



Figure 2: Convergence towards the steady state at x* =0 for: (a) 7 = 1 and; (b) 
7 = 2. The initial conditions are characterized in the figure. For this figure, we 
used the parameter = 0. 

labelled in the figure. 

We see from Fig. [2] that depending on the initial condition 
Xo, the orbit stays confined in a plateau of constant x and, af¬ 
ter reaching a crossover iteration number, Ux, the orbit suffers 
a changeover from a constant regime to a power law decay 
marked by a critical exponent p. The length of the plateau also 
depends on the initial xq. Based on the behaviour observed from 
Fig. [ 2 ] we can suppose that: 

1. For a sufficient short n, say n ^ Ux, the behaviour of x n 
is given by 

x(n) oc Xq , for n Hx , ( 2 ) 

and because x oc xq, we conclude that the critical exponent 
a = 1. 

2. For sufficient large n, i.e., n » Ux, the dynamical variable 
is described as 

x(n) oc , for n ::$> Ux , (3) 

where the exponent p is called a decay exponent. The nu¬ 
merical value is not universal and depends on the nonlin¬ 
earity of the mapping. 
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Figure 3: Plot of the crossover iteration number rix against the initial condition 
XQ together with their power law fitting for 7 = 1 with slope z = -1.0002(3) 
and 7 = 2 with slope z = - 2 . 001 ( 2 ). 

3. Finally, the crossover iteration number is given by 

oc 4 ’ 

where z is a changeover exponent. 

The exponents p and z can be obtained by considering spe¬ 
cific plots. After the constant plateau, a power law fitting 
furnishes p. Indeed for 7 = 1 (logistic map) we found 
P = -0.99981(3) while for 7 = 2 (cubic map) we obtained 
P = -0.49969(5). To obtain the exponent z we must have the 
behaviour of Ux vi-. vq, where Ux is obtained as the crossing of 
the constant plateau by the power law decay, as shown in Fig. 

0 

The slope obtained for 7 = 1, as shown in Fig. [3] is 
z = -1.0002(3) while for 7 = 2 the exponent obtained is 
z = - 2 . 001 ( 2 ). 

The behaviour shown in Fig. |2] together with the three scal¬ 
ing hypotheses allow us to describe the behaviour of v as an 
homogeneous function of the variables n and vq, when yu = 0 , 
of the type 

v(vo, n) = lx(Pxo, fin ), (5) 

where I is a scaling factor, a and b are characteristic exponents. 
Because I is a scaling factor we choose I^xq = 1, leading to 
/ = Substituting this expression in Eq. © we obtain 

v(vo, n) = x~'^^^x{\,x~^^^n ). ( 6 ) 

Assuming the term x{\,x~^^^n) is constant for n Ux and com¬ 
paring Eq. (O with the first scaling hypothesis we conclude that 
a = -1/a. Moving on and choosing fin = 1, which leads to 
I = and substituting in Eq. we obtain 

x(xo,n) = n~^^^x(n~^^^xo,l). (7) 

Again we suppose the term x(n~^^^xo, 1) is constant for n » 
nx. Comparing then with second scaling hypothesis we end 
up with p = -\/h. Finally we compare the two expressions 
obtained for the scaling factor. It indeed leads to nx = x^^. 



Figure 4: Overlap of all curves shown in Fig. |2onto a single and universal plot, 
after a convenient rescale of the axis, for both: (a) 7 = 1 and (b) 7 = 2 . 

A comparison with third scaling hypothesis allow us to obtain a 
relation between the three critical exponents a, p and z therefore 
converging to the following scaling law 


The knowledge of any two exponents allow to find the third 
one by using Eq. Moreover the exponents can also be used 
to rescale the variables x and ^ in a convenient way such that 
X ^ x/Xq and n n/x^^ and overlap all curves of x n onto 
a single and hence universal curve, as shown in Fig. lU 

Before moving to the next section and consider the dynamics 
in a more analytical way, let us discuss here the dynamics for 
ja ^ 0. This characterises the neighbouring of the bifurcation. 
The convergence to the steafiy state is marked by an exponential 
law of the type (see Refs. ll22.l^ ) 

x(n,fi) cc e-"'^ , (9) 

where r is the relaxation time described by 

T oc / , (10) 

where 5 is a relaxation exponent. Figure|5]shows the behaviour 
of T yu for two different values of 7 . 
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Grouping the terms properly we obtain the following differ¬ 
ential equation 



(13) 


Indeed the initial condition xq is defined for ^ = 0. Of course 
for a generic n we have x(n). Using these terms as limit of the 
integrals we end up with 



dx 



(14) 


After integrating Eq. (O and organising the terms properly we 
obtain the following expression 


x(n) = 


■^0 

\xlyn+\^'^ 


(15) 


Let us now discuss the implications of Eq. (fTSl) for specific 
ranges of n. We start with the case x^^yn 1, which is equiv¬ 
alent to the previous section of ^ Ux. Eor such a case we 
obtain that x{n) = vq. A quick comparison with first scaling hy¬ 
pothesis allow us to conclude that the critical exponent a = 
Second we consider the situation x^^yn » 1, corresponding to 
^ in the previous section. Eor such case we obtain that 


x{n) ^ n . 


(16) 


Comparing then this result with scaling hypothesis two of 
the previous section we conclude p = -1 /y. The last case is 
obtained when xj^yn = 1, which is the case of ^ = Ux. Then we 
obtain 


Figure 5: Plot of the relaxation to the fixed point as a function of ji in the 
logistic-like map for the exponents: (a) y = 1 and; (b) y = 2. 

A power law fitting furnishes the exponent 6 = -\ and is 
independent on the value of the parameter y. In the next sec¬ 
tion we describe how to obtain the exponents discussed in this 
section using an analytical approach. 


3. An analytical description to the equilibrium 


Let us now discuss a different approach to reach the equilib¬ 
rium. We start first with case (i), i.e., at the bifurcation point 
R = The equation of the mapping is then written as 

Xn+\ = Xn - ■ (11) 


Very near the fixed point, we suppose the dynamical variable 
X can be considered as a continuous variable. Therefore Eq. 
(HB is rewritten in a convenient way as (see also Ref. iH for 
a recent application in a 2-D mapping) 


^n +1 


^n+\ 

(n-\- 1) - n 


dx 

dn 


vT+l 


( 12 ) 


nx = xj . (17) 

A comparison with third scaling hypothesis gives us that 
^ = -y. With this procedure we obtained all the three critical 
exponents discussed in the previous section as function of the 
parameter of the nonlinearity y. Numerical simulations were 
made for several different values of y considering either odd, 
even, irrational and other set of numbers. The numerical find¬ 
ings confirm the validity of both the scaling law as well as the 
analytical procedure. 

The last point to discuss is the case ofR < Rc, i.e., imme¬ 
diately before the bifurcation. Eor this case we can rewrite the 
mapping as 

Xn+l -Xn = X„{R - 1) - , 

^n+\ - Xn dx 

{n+ 1) - n dn' 

= x{R-\)-Rx>'^'^ . (18) 


We have to emphasise that near the steady state x = 0 and 
considering y > 1, the term x^^^ goes faster to zero as compared 
with X. Then the last term of Eq. (fTSl) can be disregarded. With 
this approach we obtain the following differential equation 


dx 

dn 


= -xji , 


(19) 
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Figure 6 : (a) Convergence to the fixed point x* = 2/3 for 7 = 1 considering 
R = 3. The slope of decay obtained was ^ = -0.49939(7). (b) Overlap of the 
curves shown in (a) onto a single and hence universal plot. 




Figure 7: (a) Convergence to the fixed point = 0.4399601688 ... for 7 = 1 
considering R = I + ^^6. The slope of decay obtained was [3 = -0.498(1). (b) 
Overlap of the curves shown in (a) onto a single and hence universal plot. 


where ji = \ - R. Considering again that for ^ = 0 the initial 
condition is xq, we have to integrate the following equation 



After integration and grouping the terms we obtain 


x{n) = xoe-^^ . ( 21 ) 

Comparing this result with Eqs. ^ and (fTOl) we conclude 
that the exponent S = -1. This finding is in good agreement 
with the simulations shown in Fig. [Jl 


4. Discussions 


Convergence to the steady state at a period-doubling bifur¬ 
cation was also observed to obey a homogeneous function. In¬ 
deed to apply the formalism we have to look at the distance 
to the fixed point. Such observable can be defined as (follow¬ 
ing Refs, [ 2 ^) y(n) = F^(F^(x)) - x* where F stands for 
the mapping, m = 2 ,4, 6 ,... and v* is indeed the expression of 
the fixed point. In the previous sections, the fixed point was 


located at x* = 0 , hence it was convenient to rescale the ob¬ 
servables in terms of xq, which was the initial distance from the 
fixed point. In the period-doubling bifurcation considered here, 
X* is not zero anymore, hence we represent the distance from 
the fixed point as e. Furthermore, the scaling is dependent on e. 
Figure [ 6 | shows the convergence to the fixed point for 7 = 1 and 
R = 3. The slope obtained for the decay is jS = -0.49939(7). 
The crossover exponent was obtained as z = -2.001(4), and 
the same scaling law as obtained previously is applied here too. 
The exponents obtained for the convergence towards the fixed 
point in the period-doubling bifurcation show to be indepen¬ 
dent of the nonlinearity 7 . The scaling (see the overlap of the 
curves shown in Fig. Oa) onto a single and universal plot as 
shown in Fig. Ob)) is better seen after the transformations: (i) 
(x - X*) ^ (x - x*)/ 6 ^ where e stands for the distance of the 
initial condition to the fixed point x* and; (ii) n nle^. For 
y = 1, the first period-doubling bifurcation happens at R = 3, 
and the fixed point is x* = 2/3. For the second period-doubling, 
R = \ and the two fixed points are 


■^1,2 2 


1 + - + - V/?2 - 2 /? - 3 

R R 


( 22 ) 
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For R = 1 V 6 , the two fixed points assume the values = 

0.8499377796... and x* = 0.4399601688 .... In this cases, the 
critical exponents sltq a = 1, z = -2 and {3 = -0.498(1). 

The numerical results obtained for the parameter /? = 1 + V6 
are shown in Fig. [71 The convergence to the fixed point is 
plotted in Fig. [TJa) for different values of initial conditions, as 
shown in the figure while the overlap onto a single and universal 
plot is shown in Fig. [3b). The slope of the decay obtained is 
the same one as obtained for the first period-doubling bifurca¬ 
tion observed at 7? = 3. This result confirms the independence 
of the critical exponents a, JS and z, on the period-doubling bi¬ 
furcation, as a function of y. 

5. Conclusions 

To summarize, we have considered the convergence to the 
steady state in a family of logistic-like mappings in a bifurca¬ 
tion point near a transcritical, pitchfork, period-doubling bifur¬ 
cation, and around its neighbouring. At the bifurcation point we 
used a phenomenological description to prove that decay to the 
fixed point is described by using a homogeneous function with 
three critical exponents. The three critical exponents are related 
between themselves via a scaling law of the type z = czip. Near 
the bifurcation point the convergence to the fixed point is given 
by an exponential decay and the relaxation time is described 
by a power law of the type r oc [J. We found 6 = -\ and is 
independent on the nonlinearity of the mapping. The results 
obtained by the phenomenological approach is confirmed by an 
analytical description and is valid for any y > 1. The expo¬ 
nents obtained for the period doubling bifurcation were a = 

P = -0.5, z = -2 and are independent of the parameter y. 
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